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Atmospheric neutrinos as probes of neutrino-matter interactions 
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Neutrino oscillation experiments provide a unique tool for probing neutrino-matter interactions, 
especially those involving the tau neutrino. We describe the sensitivity of the present atmospheric 
neutrino data to these interactions in the framework of a three-flavor analysis. Compared to the 
two-flavor — Vr analyses, we find qualitatively new features, in particular, that large non-standard 
interactions, comparable in strength to those in the Standard Model, can be consistent with the 
data. The existence of such interactions could imply a smaller value of the neutrino mixing angle 
and larger value of the mass-squared splitting than in the case of standard interactions only. This 
and other effects of non-standard interactions may be tested in the next several years by MINOS, 
KamLAND and solar neutrino experiments. 

PACS numbers; 13.15.-|-g, 14.60.Lm, 14.60.Pq 



I. INTRODUCTION 

Historically, the studies of the neutrino properties 
have lead to several important breakthroughs in particle 
physics. The discovery of the neutral current processes 
played a key role in establishing the Standard Model 
(SM), while the recent confirmation of neutrino oscilla- 
tions provided the first real sign of non-standard physics. 
To this day, however, neutrino interactions remain one 
of the least tested aspects of the SM. The importance of 
their direct experimental determination, whether it leads 
to further confirmation of the SM or to a discovery of new 
physics, can hardly be overstated. 

In this paper, we explore the possibility of measuring 
neutrino-matter interactions using the atmospheric neu- 
trino data. A key feature of these data is their sensitivity 
to the interactions of the tau neutrino. Indeed, atmo- 
spheric muon neutrinos are believed to oscillate into tau 
neutrinos; modified interactions of the latter may make 
the oscillation pattern incompatible with the data. In 
contrast, accelerator-based experiments have restricted 
mainly the interactions of the muon neutrino. 

The sensitivity of the atmospheric neutrino data to 
non-standard interactions (NSI) was studied in P, 0. 13 ■ 
It was found that, in the context of a two- flavor ^ Vt 
analysis, these data are quite sensitive to NSI. The sen- 
sitivity to possible <^ conversion, however, was 
not systematically explored beyond some specific exam- 
ples pf. Solar neutrino data constrain this channel 0, 
but at the same time leave room for non-trivial effects, 
both on the solar and atmospheric neutrinos. In the solar 
case, the v^. <-!■ interaction gives rise to a new discon- 
nected solution LMA-0 Q, |^ |^, thus introducing an 



uncertainty in the determination of the solar oscillation 
parameters. Whether a combined analysis of the solar 
and atmospheric data eliminates this uncertainty is an 
important outstanding question. 

In the following, we perform the first study of the sen- 
sitivity of the atmospheric data to NSI in the three-flavor 
framework. Our goals are: first, to describe the physics 
behind the sensitivity of the data to NSI; second, to 
present the bounds on the interaction parameters, based 
on a detailed fit to the data; and finally, to discuss how 
the possible presence of modified interactions can affect 
the measurement of the oscillation parameters. We in- 
tend to publish the results of a complete scan of the pa- 
rameter space elsewhere Q. 



II. SETUP 

The atmospheric neutrino flux is comprised of both 
neutrinos and antineutrinos, with the energy spectrum 
spanning roughly five orders of magnitude, from 10~^ to 
lO'* GeV 13 . Initially created as the electron and muon 
flavor states, these neutrinos (depending on the distance 
traveled) undergo flavor oscillations. The oscillations are 
governed by the Hamiltonian H , which has a kinetic (vac- 
uum) part and an interaction part. The former is given, 
at the leading order, by 
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where we omit the smaller, "solar" mass splitting, Attiq 
(3, and the remaining two mixing angles, d\2 and ^13. In 
Eq. iQl, A = Am^/(4i?i,), Ey is the neutrino energy and 
Am^ is the largest mass square difference in the neutrino 
spectrum; the flavor basis, v^,v^j_,Vt, is assumed. The 
interaction part (up to an irrelevant constant) will be 
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parameterized as 
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Here Gf is the Fermi constant, rie is the number density 
of electrons in the medium. 

The physical origin of the epsilon contributions in -ffmat 
can be the exchange of a new heavy vector or scalar 22] 
particle. We parameterize the resulting NSI with the 
effective low-energy four-fermion Lagrangian 



-NSI 



+ h.c. (3) 

Here e^^^ i^Lff^) denotes the strength of the NSI be- 
tween the neutrinos of flavors a and P and the left- 
handed (right-handed) components of the fermions / and 
/. The epsilons in Eq. (jSJ are the sum of the contri- 
butions from electrons (e"^), up quarks (e"), and down 
quarks (e'') in matter: eap = J2f=u,d,e^Lf3^f /"^e- In 

Notice that 



JfR , 



turn, 



e{o = e{n -I- el^ and el^ = el^^ 



the matter effects are sensitive only to the interactions 
that preserve the flavor of the background fermion / (re- 
quired by coherence ^3) and, furthermore, only to the 
vector part of that interaction. 

Neutrino scattering tests, like those of NuTeV and 
CHARM 12], mainly constrain the NSI couplings of the 



muon neutrino, e.g., ] ] < 10 
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The limits they place on e^e , Cer , and err are rather loose. 
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< 3, -0.4 < 
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< 0.7, Je™] < 0.5, 
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e:;^! < 0.5 [13] • Stronger constraints exist on the cor- 
responding interactions involving the charged leptons. 
Those, however, cannot, in general, be extended to the 
neutrinos, for example when the underlying operators 
contain the Higgs fields and hence will not be con- 
sidered here. 

Given the above bounds we will set ee^i and e^^ to zero 
in our analysis. Furthermore, for simplicity we will also 
set e^r to zero. The earlier analyses of the atmospheric 
neutrino data have indicated that this parameter is 
quite constrained {efj.T< 10~^ — 10~^). Corrections due 
to non- vanishing e^,- will be described in 0]. Here, we 
have a three-dimensional NSI parameter space, spanned 
by Cee, Cer, and err- 



III. CONVERSION EFFECTS AND 
SENSITIVITY TO THE NSI 

The physics of the sensitivity of the atmospheric neu- 
trino data to eee, fer, and Ctt can be understood as fol- 
lows. The data are known to be very well fit by large- 
amplitude oscillations between the i'^ and i/r states. This 
holds both at high energy {E^, > 10 GeV), where only 
the muon neutrino flux is measured, and at lower ener- 
gies, where both the muon and electron neutrino data 



are available. These oscillations are driven by the off- 
diagonal — Vr mixing in Eq. J^l and the introduction 
of sufficiently large NSI for the tau neutrino will, in gen- 
eral, suppress that mixing. Since the vacuum Hamilto- 
nian scales as i?j7^, this suppression should be especially 
strong at high energy, in the through-going muon sample. 

As a simple illustration, consider the case when only 
trT is nonzero. Clearly, e^-r introduces a diagonal split- 
ting between the and Vr states, thereby decreasing 
the effective mixing angle in matter. The correspond- 
ing bound can be estimated by comparing the mat- 
ter term \pl trr GpUe to the vacuum oscillation term 
/S.m? / {2E^). For neutrinos going through the center 
of the Earth, the highest energy at which an oscilla- 
tion minimum occurs in the standard case is around 
i?o 20 — 30 GeV. If the matter term is sufficiently 
large, \/2 erT Gprie ^ ATO^/(2i?o), the mixing in mat- 
ter and hence the oscillation amplitude are expected to 
be suppressed. Substituting numerical values, we find a 
bound Ett^ 0.2. 

Next, we generalize this argument to the case of non- 
vanishing Cee, Eer- The matter part of the Hamiltonian 
-ffmat can be diagonalized by rotating in the Ve — Vr sub- 
space. In the new basis (i^e', i'^, i^t')j ^mat has the form 

(iia(7(Ae' , 0, Ai-')j with Ae'.r' — V^GFne{l+ fee + Err 

±^/{l+ eee - frrY + 4] Eer l^)/2. It Straightforwardly 
follows that if lAe',r'l ^ ^rr? I(2Eq), the oscillations of 
the muon neutrinos proceed unimpeded, while in oppo- 
site case, lAe',r'| ^ A"i^/(2£^o)j they are suppressed. 

It is very important to consider the intermediate 
regime, when the spectrum has the hierarchy (a) jA,-'] < 
AmV(2£;o) < JAe'l or (b) JAe'] < AmV(2i;o) < JA^-]- 
In both cases, the oscillations between j/^ and the corre- 
sponding light eigenstate are allowed to proceed while 
those between j/^ and the heavy eigenstates are sup- 
pressed. Remarkably, the resulting oscillation pattern is 
indistinguishable from the standard case at high energy, 
where only muon neutrinos are detected. 

From now on we specialize to hierarchy (a), which is 
smoothly connected to the origin eee=eeT=iTT— ((b) 
is realized only if eee + Ett is a large negative number). 
When it is satisfied, muon neutrinos oscillate into the 
state 



Vr' = —Spe"^^^ Ve +Cp V. 



(4) 



where = cos/5, Sfj — sin/?, 2^ — Arg{eer), tan2/3 

2] Cer 1/(1+ Cee - Ctt)- 

The condition jA,-'! < Am'^ /{2Eq) implies 



11+ eee + err ^ ^/{l+^^^^^^rF+M^^ \ ^ 0.4. (5) 

This equation gives our analytical prediction for the 
bound on eee,eeT,eTr- When eee=eeT= 0, it reduces to 
the bound err^ 0.2 given above. 

The region Eq. (|SJ) describes extends to large values of 
eer, err- To SBC this, notc that in the limit At' = 0, or 



err= 1 Eer 1^/(1+ Cee) 



(6) 
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the Vfi-^ Vt-i oscillations, though dependent on the matter 
angle /3, are independent of the absolute size of the NSI. 
As already mentioned, these oscillations have the same 
dependence on the neutrino energy and on the distance 
L as vacuum oscillations and therefore mimic their effect 
in the distortion of the neutrino energy spectrum and of 
the zenith angle distribution. More specifically, we get 
the oscillation probability: 

Pivy.^ vr') = sin^ 20™ sin2[Am^L/(4£;,)] , (7) 

where the effective mixing and mass square splitting are 
derived to be 

Am^„ = Am^ [(c2.(l + - s^f/A + [s^ecpf] , 
tan26l,„ = 2s2e 0^3/(029(1 + cj) - s}) . (8) 

If NSI are present, but not included in the data analy- 
sis, a fit of the highest energy atmospheric data, i.e. the 
through-going muon ones, would give Am^ and 9„i in- 
stead of the corresponding vacuum quantities. If we fix 
a set of NSI and - to reproduce the no-NSI case - re- 
quire that 9m — 7r/4 and Am^ ~ 2.5 • 10~^ eV^, from 
Eqs. ||SJ) we get that the vacuum mixing would not be 
maximal; in particular we have cos 20 ~ s^/(l -I- c^) and 
Am^ ~ Am^(l + cos"^ /3)/2. 

In the intermediate energy range, E ^ 1 — 10 GeV, 
when matter and vacuum terms are comparable, the re- 
duction to a two-neutrino system is not possible, and the 
problem does not allow a simple analytical treatment. 
The neutrino conversion probability in this energy range 
depends on the sign of the neutrino mass hierarchy (nor- 
mal, Ato^ > 0, or inverted, Ato^ < 0). At the sub-GeV 
energies, we expect vacuum-domination, and therefore 
small deviations with respect to vacuum oscillations . 

Finally, we observe that for 6*13 = 0, Atoq — 0, as has 
been assumed here, there is no sensitivity to ip, the phase 
of 1^ • This is unlike the case of the solar neutrinos, 
where -0 plays a crucial role 0. Corrections due to ^13 
and Atuq 7^ break the phase degeneracy and will be 
presented elsewhere 0- 



IV. NUMERICAL RESULTS 

We performed a quantitative analysis of the atmo- 
spheric neutrino data with five parameters: two "vac- 
uum" ones, (Am^,0), and three NSI quantities (eee,eer 
,err)- The goodness-of-fit for a given point is deter- 
mined by performing a fit to the data. We use the 
complete 1489-day charged current Super-Kamiokande 
phase I data set [13 , including the e-like and /i-like data 
samples of sub- and multi-GcV contained events (each 
grouped into 10 bins in zenith angle) as well as the stop- 
ping (5 angular bins) and through-going (10 angular bins) 
upgoing muon data events. This amounts to a total of 
55 data points. For the calculation of the expected rates 
we use the new three-dimensional atmospheric neutrino 



fluxes given in Ref. The statistical analysis of the 

data follows the appendix of Ref. Q . 

The results of the K2K experiment have been included. 
Their addition has a minimal impact on our results, pro- 
viding some constraint at high Am^ . The details of the 
K2K analysis can be found in Sec. 2.2 of Ref. [T^ . 
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FIG. 1: A 2-D section (eee= —0.15) of the allowed region of 
the NSI parameters (shaded). We assumed Attiq = and 
613 = 0, and marginalized over 6 and Am^. The dashed con- 
tours indicate our analytical predictions. See text for details. 

Upon scanning the parameter space and marginalizing 
over Ato^ and 9 we obtain the three-dimensional allowed 
region in the space (cee, Cer, Crr)- As an illustration, in 
Fig. n we show a section of this region by the plane eee= 
—0.15 (the choice motivated by the solar analysis in jj]). 
The minimum occurs at eeT= 0.07, 6,-,-= 0.01; the 
value at the minimum, Xmin — 48.50, is virtually the 
same as at the origin (no NSI), Xorig = 48.57. The shaded 
regions correspond, from the innermost contour, to ~ 
Xmin < 7.81, 11.35, and 18.80. They represent the 95%, 
99%, and 3.6(t confidence levels (C.L.) for three degrees 
of freedom (d.o.f.). The last contour also corresponds to 
the 95% C.L. for 50 d.o.f.. For the purpose of hypothesis 
testing this means that a theory which gives NSI outside 
of this region should be rejected. 

The dashed-dotted parabola illustrates the condition 
of zero eigenvalue, Eq. ©; the two outer curves give 
the predicted bound according to Eq. jSJ). For both, 
the agreement between the theory and numerical results 
is quite convincing. Moreover, we have verified that 
the agreement remains very good for e^e in the range 
—0.7 <eee< 0.3 For the case when only Cer is non- 
zero we find the bounds | Cer \ < 0.38 at 99% C.L. and 
I Eer I < 0.5 at 3.6a. 

The extent of the allowed region along the parabola is 
beyond the scope of our analytical treatment. Indeed, 
since at high energy the leading NSI effect is canceled by 
construction, the fit quality is determined by subdomi- 
nant NSI effects in all energy samples. Remarkably, these 
effects are rather small, especially for the inverted mass 
hierarchy, where the region x^ — Xmin — "^-^^ extends up 
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FIG. 2: The effect of the NSI on the aUowed region and best- 
fit values of the oscillation parameters; see text for details. 



to eer= li as the Figure shows. The extent is somewhat 
smaller (to eer= ±0.6) for the normal hierarchy 0- 

As mentioned, the symmetry of Fig. ^ with respect to 
follows from setting ^13 and Am^ to zero. 
We have checked that for 6*13 ^ the region does indeed 
develop a small asymmetry 7]. 

Finally, we notice that as the NSI point is moved 
along the parabola © away from the origin, the best- 
fit oscillation parameters change. This behavior is illus- 
trated in Fig. El The four contours plotted are the sec- 
tions of the five-dimensional allowed region in the space 
(Am^, 9, eee, egr, Ett) for the values: (i) £ee— —0.15, eer— 
0,err^ 0; (ii) eee= -0.15, eer= 0.30, e^^= 0.106; (iii) 
eee= -0.15, eer= 0.60, 6^^= 0.424; (iv) eee= -0.15, eer = 
0.90, eTT= 0.953. The contours are drawn for Ax^ = 
11.07 with respect to the global minimum, marked by 
a star (this corresponds to 95% CL for 5 d.o.f.). The 
dots represent the best fit points obtained in each of the 
four cases above, when the NSI are assumed known. The 
figure confirms the trend predicted in Eq. (|SJ|: as the 
NSI increase the best-fit mixing angle becomes smaller, 
9 < tt/4, while the best-fit Am^ increases. 

The shaded region in the background represents the 
result of marginalizing with respect to eer and err (for 
eee= —0.15). It was drawn for Ax^ = 11.83, correspond- 
ing to 99.73% (3cr) CL. for 2 d.o.f. While the best-fit 
point is nominally unchanged by the NSI, the shape of 
the function is significantly altered, so that small val- 
ues of the mixing angle are now much less disfavored. 

V. CONCLUSION AND IMPLICATIONS 

In conclusion, the three-flavor oscillation analysis ex- 
hibits qualitatively new features compared with the <-> 
Vr studies. While giving no positive evidence for non- 
standard physics, the data do not exclude large NSI. 
The old bound \ e^r \ < 0.15 at 99% CL. is signifi- 
cantly relaxed by the addition of nonzero eer- if fer and 
Err are varied simultaneously, even erT~ 1 is allowed. 



The effect is even more dramatic when one compares the 
bounds on the flavor-changing interactions: the bound 
ep_T< a few x 10~^ Q does not extend to e^r, for which 
order one values are allowed. The bound eer< 0.38 at 
99% CL. - found for the case when the other epsilons 
vanish - is likewise new. We also stress that the effects 
of these large NSI are absolutely non-trivial. In particu- 
lar, the best-fit values of the oscillation parameters shift 
to smaller angles and higher Am^. The physical mecha- 
nism behind this is well understood analytically. 

These properties have important implications for other 
neutrino experiments. For example, the MINOS exper- 
iment will measure the mixing angle 9 with good pre- 
cision; assuming the true value 9 ^ tt/A one expects 
sm^9 > 0.4 at 99% CL. [l|. Given its baseline of 
only 735 km, MINOS will be almost free from matter 
effects and will therefore determine the vacuum value of 
9. Thus, it can look for the large NSI effects described 
here: a measurement of sin^ 9 > 0.4 and low Am^ would 
contribute to disfavor large NSI, while in the opposite 
case (sin^e ~ 0.3) one would have indications of non- 
standard physics. 

Further tests can be performed at future experiments 
with longer baseline (a few x lO'^ km) and sufficiently high 
energy, E^, ~ 20—30 GeV. These will be directly sensitive 
to the NSI matter effects and a discrepancy between them 
and MINOS could be a sign of NSI. 

It is worth mentioning that, in principle, another sig- 
nature of large NSI could be an unexpected number of 
neutral current events in a detector relative to the num- 
ber of charged current events for a given oscillation sce- 
nario. This effect would be due to the modification of the 
detection neutral current cross section by large eer and 
Err- Unfortunately, no firm prediction on the size of the 
effect can be made, since the effect additionally depends 
on the values of the axial coup lings that do not enter the 
propagation Hamiltonian |25l |. 

Finally, there are important implications for solar neu- 
trinos. Indeed, Cee, Cer, and Crr also enter the solar neu- 
trino survival probability, and may alter the interpreta- 
tion of the current data 4,^,i6||. Our results here demon- 
strate that the LMA-0 solution, obtained with non-zero 
Eer as in 01, is compatible with the atmospheric data. 
On the other hand, solar neutrinos, being sensitive to 
the phase of eer, rnay rule out a sizable portion of the 
parameter space with eeT< 0, 0- This exclusion is 
impacted by the latest results from KamLAND j3| and 
it would be very interesting to see the combined bound 
from the atmospheric, solar, and KamLAND data. 



Acknowledgments 

We are very grateful to Thomas Schwetz for sharing 
with us his K2K analysis. We thank the Aspen Center 
for Physics where part of this research was completed 
for hospitality. A. F. was supported by the Department 
of Energy, under contract W-7405-ENG-36, C. L. by a 



5 



grant-in-aid from the W. M. Keck Foundation and the PHY-0354776. 
NSF grant PHY-0070928, and M. M. by the NSF grant 



[1] N. Fornengo, M. Maltoni, R. T. Bayo and J. W. F. Valle, 
Phys. Rev. D 65, 013010 (2002) arXiv:hep-ph/0108043 . 

[2] M. Guzzo, P. C. de Holanda, M. Maltoni, H. Nunokawa, 
M. A. Tortola and J. W. F. Valle, Nucl. Phys. B 629, 
479 (2002) arXiv:hep-ph/0112310 . 

[3] M. C. Gonzalez-Garcia and M. Maltoni, 
arXiv:liep-ph/0404085 (to appear in PRD). 

[4] A. Friedland, C. Lunardini and C. Pena-Garay, Phys. 
Lett. B 594, 347 (2004) arXiv:hep-ph/0402266 . 

[5] M. M. Guzzo, P. C. de Holanda and O. L. G. Peres, Phys. 
Lett. B 591, 1 (2004) arXiv:hep-ph/0403134 . 

[6] O. G. Miranda, M. A. Tortola and J. W. F. Valle, 
arXiv:hep-ph/0406280 

[7] A. Friedland, C. Lunardini, M. Maltoni, in preparation. 

[8] R. Engel, T. K. Gaisser and T. Stanev, Phys. Lett. B 
472, 113 (2000) arXiv:hep-ph/9911394 . 

[9] C. W. Kim and U. W. Lee, Phys. Lett. B 444, 204 
(1998); O. L. G. Peres and A. Y. Smirnov, Phys. Lett. 
B 456, 204 (1999); M. C. Gonzalez-Garcia, M. Maltoni 
and A. Y. Smirnov, arXiv:hep-ph/0408170 
[10] A. Friedland and C. Lunardini, Phys. Rev. D 68, 013007 

(2003); JHEP 0310, 043 (2003). 
[11] G. P. Zeller et al, Phys. Rev. Lett. 88, 091802 (2002) 

[Erratum-ibid. 90, 239902 (2003)] hep-ex/0110059 . 
[12] P. Vilain et al., Phys. Lett. B 335, 246 (1994). 
[13] S. Davidson, C. Pena-Garay, N. Rius and A. Santamaria, 

JHEP 0303, Oil (2003) arXiv:hep-ph/0302093 . 
[14] Z. Berezhiani, A. Rossi, Phys. Lett. B 535, 207 (2002). 
[15] Y. Hayato, Super-Kamiokande Coll., in Pro- 
ceedings of the International Europhysics Con- 
ference on High Energy Phystcs, Aachen, July 
2003, [Eur. Phys. J. C 33, s829 (2004), supple- 
ment] , http : / / eps2003 . physik . rwth-aachen . de/ 



transparencies/07/ index .php. 

[16] M. Honda, T. Kajita, K. Kasahara and S. Midorikawa, 
astro-ph/0404457 

[17] M. Maltoni, T. Schwetz, M. A. Tortola and 
J. W. F. Valle, arXiv:hep-ph/0405172 

[18] V. D. Barger, A. M. Gago, D. Marfatia, W. J. C. Teves, 
B. P. Wood and R. Zukanovich Funchal, Phys. Rev. D 
65, 053016 (2002) arXiv:hep-ph/0110393j. 

[19] T. Araki et al., arXiv:hep-ex/0406035 

[20] S. Fukuda et al., Phys. Rev. Lett. 85, 3999 
(2000) arXiv:hep-ex/0009001 ; H. Sobel, in proceed- 
ings of Neutrino 2000, XlXth International Confer- 
ence on Neutrino Physics and Astrophysics, http:// 
nu2000 . sno . laurentian. ca/H. Sobel/ index .html; M. 
Shiozawa, in proceedings of Neutrino 2002, XXth 
International Conference on Neutrino Physics and 
Astrophysics, http : //neutrino2002 .ph. turn, de/pages/ 
transparencies/ shiozawa/ index . html. 

[21] S. N. Ahmed et al. [SNO Collaboration], Phys. Rev. Lett. 
92, 181301 (2004) arXiv:nucl-ex/0309004 . 

[22] Scalar interactions of the form (f'./fl )(,/«!') reduce to 
Eq. (|3J upon the application of the Fierz transformation, 
see e.g. S. Bergmann, Y. Grossman and E. Nardi, Phys. 
Rev. D 60, 093008 (1999) arXiv:hep-ph/9903517 . 

[23] The matter effects may give beat-like structure to the 
survival probability as a function of energy E„, which 
tends to become unobservable upon integration over E^. 

[24] Clearly, by redefining the phases of and Ut one changes 
the phases of the ejj. and er entries of the Hamiltonian. 

[25] It is for this reason that the NC information .20. . .2L1 was 
not used in our analysis here and in 4]. 



